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Abstract. A new approach to the analysis of three body decays is presented. Measurements of the 
5-wave Kn amplitude are made in independent ranges of invariant mass from threshold up to the 
upper kinematic limit in D + — ► K~n + K + decays. These are compared with results obtained from 
a fit where the 5-wave is assumed to have k and ^(1430) resonances. Results are also compared 
with measurements of K K + elastic scattering. Contributions from l= \ and l= \ are not resolved 
in this study. If / = j dominates, however, the Watson theorem prediction, that the phase behavour 
, below Kn' threshold should match that in elastic scattering, is not well supported by these data, 

■"sj- ■ Production of K k + from these D decays is also studied. 

o 
o 

in : INTRODUCTION 
O 

><! \ Decays of D and B mesons show promise as a source of information on the light-quark 

V ' mesons they produce. Their decays to S- wave systems in three pseudo-scalar final states, 

^h. may help to improve our knowledge of the particularly confusing scalar meson (J p = + ) 

spectrum. Until now, extracting such information has been done in model-dependent 
ways that make assumptions about the scalar states observed that can influence the 
results. With large, clean samples of such decays, anticipated to be coming from the 
B factories and the Tevatron collider, the need for new approaches is a priority. 

Knowledge of strange scalar mesons has relied on measurements of 5-wave K n + 
scattering. These come principally from SLAC experiment El 35 (LASS) [1] and cover 
the invariant mass range above 825 MeV/c 2 . Data from other experiments below this 
range exist, but with less precision [2]. More measurements in the low mass region are 
required if the possibility of the existence of a K state is to be properly evaluated. 

In this paper, a study of the decays [Note] D + — > K~k + k + observed in data from 
Fermilab experiment E791 is presented. In an earlier analysis [3] the K~n + S- wave was 
modelled with Breit-Wigner (BW) amplitudes for K and ^(1430) resonances. A less 
model-dependent analysis is presented here. 5-wave amplitudes are measured without 
assuming any specific dependence on K n + invariant mass, or on the presence of scalar 
resonances. Results are compared with the amplitude from our earlier analysis, and also 
with measurements from LASS. 
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DATA SAMPLE 



The selection process for events used in this paper is described in Ref. [3]. A signal 
consisting of 15,079 D + — > K~7l+n£ decays, with a purity of ~ 94%, is obtained. 
Fig. 1 shows the Dalitz plot with the K~n+ invariant mass squared plotted vs. that 
for the K~Ttu system. Horizontal (and symmetrized vertical) bands corresponding to 
the A* (892) resonance are clearly seen. A complex pattern of both constructive and de- 
structive interference is seen near 2 (GeV/c 2 ) 2 due to presence of £Jf(1430), ^(1410) 
and A"| ( 1 430) . Evidence for Kl ( 1 680) , difficult to see due to smearing of the Dalitz plot 
boundary resulting from the finite resolution in the three-body D + mass, may also exist. 
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FIGURE 1. Dalitz plot forZ) + — > K~n][n£ decays. The squared invariant mass of one K~K + combina- 
tion is plotted against the other combination. The plot is symmetrized, each event appearing twice. Lines 
in both directions indicate values, equally spaced in squared effective mass, at each of which the 5-wave 
amplitude is determined by the method described in the text. Kinematic boundaries for the Dalitz plot are 
drawn for three-body mass values M = 1.810 and M = 1.890 GeV/c 2 , between which data are selected. 

An asymmetry along the K* (892) bands, most easily described by interference with 
a significant 5-wave component to these decays, is also observed. Information on this 
5-wave amplitude is obtained from its interference with the K* (892), and also the other 
well-established resonances in the Dalitz plot. 



METHOD 

In Ref. [3], as in most earlier analyses of D decays to three pseudo-scalar particles 
ijk, the isobar model, with BW resonance forms, is used. The decay amplitude si is 
described by a sum of quasi two-body terms D — > R + k, R^> i + j, in each of the three 
channels k = 1,2,3: 

si = d e lS ° + jrd n e iS » 2 Fr{ ~ P ^ J 1 — X F D (q,r D ,J)Mj(p,q) (1) 
„=i mi n -s-im Rn T Rn {s) 



The squared invariant mass of the ij system is s, 7 is the spin, m Rn the mass and T Rn (s) the 
width of each of the N resonances R n seen to be contributing to the decay. For 7 > 0, F R 
and F D are Blatt-Weisskopf form factors [4], with effective radius parameters r R and r D , 
for all R n and for the parent D meson, respectively. For 7 = 0, a Gaussian form suggested 
by Tornqvist [5] is used for the D. The momenta, p and q for i and k, respectively, are 
defined in the ij rest frame, and Mj(p,q) = (—2pq)Pj(p ■ q) is introduced to describe 
spin conservation in the decay. The complex coefficients d n e lS " (n = 0,N) depend on 
the D decay and are determined by a fit to the data. The first term describes non- 
resonant (NR) decay to 1 + j + k with no intermediate resonance, and is assumed to be 
independent of s. For D + — > K^tz^k^ decays we Bose-symmetrize srf with respect to 
interchange of it^ and %t . 

In Ref. [3], an excellent fit to the data is obtained with K and ^(1430) resonant terms 
and the NR term comprising the 7 = part of Eq. (1). 

In this paper, the K 7t + S-wave is examined in a less model-dependent way, outlined 
here. A more detailed account is given in [6]. 

Terms appearing in Eq. (1) are grouped according to the value of 7. The S- wave part 
(all terms with 7 = 0, including the NR term) is factored 

y = S(s)xM (p,q)F D (q,r D ,0) (2) 

into a partial wave S(s), describing K 7t + scattering, and the product Mo(p,q)F D (q, r D ) 
describing the D decay. The P- and D- (7 = 1,2, respectively) waves are factored in the 
same way: 

& = P(j) xMi(p,q)F D (q,r D ,l) ; 9 = D(s) xM 2 (p,q)F D (q,r D ,l), (3) 

with partial waves P(s) and D(s) consisting of resonant terms given as in Eq. (1). The 
S-wave S(s), however, is replaced by a set of values C££ i% defined at 40 discrete squared 
invariant masses s = s^. These are indicated by the lines in Fig. 1. The Ck and % values 
are regarded as independent parameters to be determined by the data. 

An unbinned likelihood fit is made to the data, in a similar way to that described 
in Ref. [3]. An incoherent function describing the 6% background in the sample from 
events that are not true D decays is added in proper proportion at each three-body 
K~iz + iz + mass to a signal distribution proportional to \s^\ 2 . There are 86 parameters - 
all (cjt, J k ) and the coefficients d k e iS * for K* (892) , K{ ( 1 680) in the P-wave and K* ( 1430) 
in the D-wave. For the AT* (892), di c e' Sk = 1 is used to provide the reference amplitude. 

RESULTS AND COMPARISON WITH MODEL-DEPENDENT FIT 

This fit also results in an excellent description of the data. Comparison of the observed 
and predicted population of the Dalitz plot gives a j} probability of 50% for 363 bins. 
The S-, P- and D-waves resulting from the fit are shown in Fig. 2. 

These results are compared with the model-dependent fit from Ref. [3] which pro- 
vided an excellent description of the observed Dalitz plot distribution. Partial waves for 
this fit are also shown in Fig. 2. The main S-wave features of both fits agree well. Some 



differences, particularly at the highest and lowest ends of the invariant mass range, re- 
sult from shifts in the P- and D- waves. Resonant fractions and the total S- wave fraction 
(about 75%) also agree well within statistical limits. 
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FIGURE 2. (a) Phases Yk = ^S^^k) an d 0>) magnitudes = |5(s^)| of 5-wave amplitudes for K~n + 
systems from D + — > K~k + k + decays with the amplitude and phase of the /T*(892) as reference. Solid 
circles, with error bars, show the values obtained from the model-independent fit described in the text. The 
effect of adding systematic uncertainties in quadrature is indicated by extensions on the error bars. The 
magnitudes plotted include a D + form-factor Fo(q,ro,0) = e~ 25( *> 12 . The phase and magnitude of P(s) 
are shown in (c) and (d), respectively. In each case, the parameters and error matrix from the fit are used 
in Eq. (3) to produce the solid curves shown one standard deviation above and below the central values. 
Similar plots for the D-wave amplitude D(s) are shown in (e) and (f). In all plots, the dashed curves show 
similar, one standard deviation limits for the amplitudes obtained from the isobar model fit in Ref. [3]. In 
(a), (c) and (e), / = A phases for K k + scattering measured in the LASS experiment are shown as x's 
with error bars indicating statistical uncertainties. 

It can be concluded that, with the present sample size, no significant distinction 
between the model-dependent and model-independent parametrizations of the 5-wave 
can be made. 

A comparison of the S- wave amplitudes S(s) measured here with the amplitudes 
T(s) obtained in K~n + scattering is now made. For each partial wave J (for each iso- 
spin I) it is expected that S(s) = Q(s)T(s)/F D (q,r D ,J))*/s/p( J+ V where Q(s) describes 
the ^-dependence of K~n + production in D decays. The Watson theorem [7] requires, 
provided there is no re-scattering of the K^K^ from iG ', that Q(s) is a real function, so 
that phases found in D decay should match those in K~n + elastic scattering data. 

1=1/2 phases measured by LASS are plotted in Fig. 2. There is a large offset in the 
5-wave, about 75°, not seen in P- or D-waves. The shapes of S- and P-waves are also 
not the same. Unless significant admixture of / = 3/2 K n + production occurs, these 



results suggest that the conditions for the Watson theorem are not met in these data. 

The production function Q(s) is shown in Fig. 3. The scattering amplitude T(s) is 
assumed to be elastic, T = sin(y— 70). The phase offset Jq is required to account for 
the difference between elastic scattering and D decay. There appears to be significant 
^-dependence above ~ 1.2 GeV/c 2 that grows in the region near 1.4 GeV/c 2 where 
7~7o- 
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FIGURE3. The quantities p/\f{sk) X |S(sjt)| X F D (q,r D , 0)/sin(%— 7o) plotted as solid circles for each 
point obtained for the 5-wave amplitude in the model-independent fit described in the text. Three points 
between 1400 and 1450 MeV/c 2 are omitted from the plot as their values for sin(/£ — 70) are very small, 
making their errors extremely large. The region between the dashed lines shows the one standard deviation 
limits of this quantity for the 5-wave amplitude obtained from the model-dependent fit. The inset shows, 
as small open circles, the quantities |7o(s)|/ sin(y) taken from the LASS experiment. The elastic range, 
up to Krj' threshold, is indicated by double-headed arrows in both plots. 
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